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We present a search for the pair production of first generation scalar leptoquarks {LQ) in data 
corresponding to an integrated luminosity of 5.4 fb~^ collected with the DO detector at the Fermilab 
Tevatron Collider in pp collisions at y/s — 1.96 TeV. In the channel LQLQ — >■ eqv^q , where q,q 
are u or d quarks, no significant excess of data over background is observed, and we set a 95% C.L. 
lower limit of 326 GeV on the leptoquark mass, assuming equal probabilities of leptoquark decays 
to eq and v^q ■ 

PACS numbers: 13.85.Rm, 14.80.Sv 



Because of the limitations of the standard model (SM), 
several extensions have been proposed, among them su- 
persymmetry (SUSY), grand unified theories, and string 
theory. Many of these extensions predict the existence 
of particles that directly connect the lepton and quark 
sectors. By combining leptons and quarks in multiplets 
of a larger symmetry group, they are expected to in- 
teract through new mediating bosons called leptoquarks 
{LQ) [l|, [3|. Leptoquarks can be either scalar or vec- 
tor fields. This Letter will focus on the search for scalar 
leptoquarks, and in the following we will not distinguish 
particles from antiparticles. Because effective models are 
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assumed for experimental searches for leptoquarks, our 
search is independent of specific extensions of the SM. 

In pp collisions such as at the Tevatron Collider, lep- 
toquarks can be produced in leptoquark-antileptoquark 
pairs. Leptoquark pair production can occur via both 
quark-antiquark annihilation and gluon-gluon fusion, al- 
though quark-antiquark annihilation is dominant. The 
production cross section for scalar leptoquarks depends 
only on the strong coupling constant and on the lep- 
toquark mass, and is known at next-to-leading order 
(NLO) Q. 

Once produced, leptoquarks can decay to two final 
states: Iq and vq' (where I = e, /i, or r). It is assumed 
that in the low energy limit there is no intergenerational 
mixing. For first generation LQ pairs the final state will 
contain a pair of leptons (e or v^) and a pair of quarks 
(u or d). In this Letter, the case in which one leptoquark 
decays to eq and the other to I'eq' is considered (charge 
conjugate states are assumed in the Letter). 

We define /3 to be the branching ratio of a first genera- 



tion leptoquark to decay to eg. Then the probabihty for 
a leptoquark to decay to Vi^q' is (1 — /3), and the prob- 
abihty for a leptoquark pair to decay to the final state 
eqveq' is BR{LQLQ -^ eqVeq') = 2/3(1 - 13). Thus, the 
probability for the final state eqveq' is maximized when 
^ = 0.5. 

Limits on the production of first generation lepto- 
quarks have been reported by the DELPHI |4j] , OPAL [5i, 
H, HI 0, ZEUS S, CDF 0, and DO [3 Collaborations. 
Recently, CMS [ill [13, and ATLAS ul published the 
first searches for scalar leptoquark pair production at the 
CERN LHC. 

The DO detector consists of tracking, calorimeter, and 
muon systems [l4l - [l6| . The central-tracking system 
consists of a silicon microstrip tracker and a central 
fiber tracker, both located within a 2 T superconducting 
solenoid. A liquid-argon and uranium calorimeter con- 
sists of a central section (pseudorapidity \ri\ < 1.1 jl7|) 
and two end sections (1.5 < |?7| < 4.2). The calorimeters 
have fine transverse and longitudinal segmentation with 
three principal layers identified as electromagnetic, and 
fine and coarse hadronic. An outer muon system (|?7| < 2) 
consists of a layer of tracking detectors and scintillation 
trigger counters in front of 1.8 T toroids, followed by two 
similar layers after the toroids [l3|. Data were collected 
with the DO detector at the Fermilab Tevatron pp Col- 
lider operating at ^/s = 1.96 TeV between August 2002 
and June 2009, and correspond to an integrated luminos- 
ity of 5.4 fb"!. 

An electron is identified from energy deposits in the 
electromagnetic calorimeter that are consistent with the 
shower development expected for an electron and have a 
matching track extrapolated from the central tracker. 

Jets are reconstructed using a midpoint cone algo- 
rithm, with a cone size of 0.5 [13] . The jet energy is 
corrected to the particle level using jet energy scale cor- 
rections determined from data [2(|. The missing trans- 
verse energy (^t) is reconstructed from all the cells of the 
electromagnetic and hadronic calorimeters, except for the 
coarse hadronic sector where a noise-reduction algorithm 
is applied. Additional corrections are then applied for all 
identified objects including jets, electrons, and muons. 

Events must satisfy at least one trigger from the single- 
electron and electron-f jets suites of triggers. For aU data 
samples, trigger objects are required to match the recon- 
structed objects. The trigger efficiencies are measured in 
data and parameterized for specific lepton and jet iden- 
tification criteria. 

Scalar leptoquark pair Monte Carlo (MC) samples are 
generated using pythia [21| with CTEQ6L1 ^^ par- 
ton density functions. Signal samples are produced for 
different LQ masses between 200 and 360 GeV. The cor- 
responding cross sections at NLO are listed in Table HI 

Diboson (TI^VF, WZ and ZZ) background samples are 
produced with pythia making use of the parton distri- 
bution functions CTEQ6L1. The ti and V{V = W or 
.Z)-|-jets events are simulated with the matrix-element 
generator ALPGEN [23I|, interfaced to pythia for subse- 



TABLE I: Scalar LQ pair production cross sections, calcu- 
lated at NLO, for different Mlq M- 



Mlq (GeV) 


200 


210 


220 


230 


240 


250 


260 


270 


a(fb) 


268 


193 


141 


103 


76 


56 


42 


31 


Mlq (GeV) 


280 


290 


300 


310 


320 


340 


360 




a(fb) 


23 


17 


13 


10 


7.4 


4.2 


2.4 





quent parton showering and hadronization. Single top 
quark production is simulated using COMPHEP [2J] . The 
cross sections for background processes are calculated 
at NLO (diboson [25[)and next-to- next-to-leading or- 
der (NNLO) (F-hjets ^] and tt [21]). We correct the 
generated spectrum of the transverse momentum {pt) of 
the Z boson in MC to match a corresponding dedicated 
measurement [2g|. The pt spectrum of the W boson 
is corrected taking into account the differences between 
predicted Z and W boson pT spectra at NNLO [2^ . 

A full GEANT-based detector simulation program [30| . 
followed by the same reconstruction program as utilized 
for data, is used to process signal and background events 
from MC. In order to model detector noise and contri- 
butions from the presence of additional pp interactions, 
events from randomly selected beam crossings with the 
same instantaneous luminosity profile as data are over- 
laid on the simulated events. Background from multijet 
production (MJ), where one of the jets mimics an elec- 
tron, is evaluated from data using a data driven tech- 
nique [3l|. In MC simulations, electron energies are cor- 
rected so they match those obtained from data. In addi- 
tion, residual differences in jet energy scale and resolution 
between data and MC are reduced by applying dedicated 
corrections to MC events. 

In the eqVf^q' final state, it is not known a priori how 
to assign the jets to the leptoquark decaying to eq or 
Veq' ■ Therefore, to reconstruct the properties such as 
mass and pT of the leptoquarks from the final products, 
an algorithm is needed to choose the best pairing. We do 
not impose a requirement on the number of jets, but we 
use only the two leading px jets for pairings. There are 
two possible combinations, corresponding to the leading 
jet pairing with either the electron or the neutrino. We 
found that it is most effective to choose the pairing that 
minimizes the difference between the transverse masses, 
Mt — \JE\ — p^, where Et and pr are the transverse 
energy and the transverse momentum vector of the two 
LQs. This pairing algorithm is successful in making the 
correct assignment in about 75% of MC signal events. 

Events are selected to be consistent with the LQLQ ^■ 
eqveq' process. We require one electron with px > 
15 GcV in the central calorimeter region |77e| < 1.1; 
fx > 15 GeV, to be consistent with the undetected 
neutrino; and at least two jets with pT > 20 GeV and 
|77jct| < 2.5. To suppress MJ background, events are re- 
quired to satisfy |:T/50-hM|,''/70 > 1, where M^" is the 
transverse mass of the (e, v) combination, and ^t and 
M^" are in GeV. 



TABLE II: Event counts and the predicted number of signal 
events for Mlq ~ 260 GeV and /3 = 0.5 after each selection 
requirement. 





Data 


Total background 


Signal 


Preselection 


65992 


65703 ± 5958 


50 ± 7 


M^"" > 110 GeV 


990 


986 ± 82 


34± 5 


J2Mlq >350 GeV 


64 


55 ± 4 


27± 4 


St > 450 GeV 


15 


15 ± 1 


24± 3 



At this stage we observe 65992 data events, while we 
expect 65703 ± 61 (stat)±5958(sys) from SM background 
and 50.4±0.4(stat) ±6.8(sys) events from scalar LQ pro- 
duction for Mlq = 260 GeV and /3 = 0.5. Figure [Ila) 
shows the M^'^ distribution for the data and SM pro- 
cesses. Data are consistent with the SM predictions. To 
reduce the dominant SM y-f-jets background, we require 
My" > 110 GeV. The pairing algorithm described previ- 
ously allows us to reconstruct Mlq. Since the longitu- 
dinal component of the neutrino momentum, pz, is not 
measurable, we reconstruct only the visible mass of the 
decay LQ — t- Veq' as Mlq = M{ict + t^vis), where the 
four vector of i^vis is given as {fix,i)y, 0, fr)- Figure [IJb) 
shows the distribution of the sum ^ Mlq of the invariant 
mass of the decay LQ — > eq and the visible mass of the 
decay LQ — >■ v^q' after the requirement M^ > 110 GeV. 
We then use ^Mlq to reduce SM backgrounds, fur- 
ther requiring that ^ Mlq > 350 GeV. Finally, we re- 
quire that the scalar sum of the pt of the lepton, the 
fx, and the two jets, St, shown in Fig. [TJc) after all 
selections, be greater than 450 GeV. Selection criteria 
are optimized to achieve the best expected sensitivity for 
Mlq = 260 GeV. This yields 15 observed events for an 
expected background of 14.8±0.6(stat)±l.l(sys) events. 
The event counts after each requirement are shown in Ta- 
ble HI 

Systematic uncertainties which affect only the nor- 
malization of the background and the signal efficiency 
include uncertainties on cross sections of signal (10%) 
and background (6% — 10%) processes, normalization of 
the MJ background (20%), integrated luminosity (6.1%), 
and lepton trigger and identification (4%). Uncertain- 
ties which also affect the differential distribution of St 
which is the quantity used to set the limits on LQ are due 
to the jet energy resolution and scale, jet identification 
efficiency, parton distribution functions, and the model- 
ing of the jet Pt distribution of the dominant W^+jets 
background. Their impacts are evaluated by repeating 
the analysis with values varied by ±1 standard deviation 
(SD). For the uncertainty on the jet pt modeling, the im- 
pact is estimated by comparing the jet pt distributions 
between ALPGEN and data unfolded to particle level from 
the recent DO measurement [32]. The ratio is applied as 
weight to the M^-|-jets jet pt distribution, and the new 
distribution is taken as ±1 SD band. 

The distribution of the St after all selection require- 
ments, shown in Fig. [IJd), is used as a discriminant to 
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FIG. 1: (color online) (a) Mf." distribution after preselec- 
tion, (b) J2Mlq for M!r" > 110 GeV, (c) the St for 
M^r" > 110 GeV and J2Mlq > 350 GeV, and (d) the St 
distribution after the final selection, which is used to set an 
upper limit on the LQ pair production cross section. 



set an upper limit on the LQ pair production cross sec- 
tion in the eqv^q' channeL For each generated Mj^q, 
the Hmit is calculated at the 95% C.L. using the semi- 
frequentist CLg method based on a Poisson log-likelihood 
test statistic [33| . Signal and background normalizations 
and shape variations due to systematic uncertainties are 
incorporated assuming Gaussian priors. The best fit to 
the background distributions is evaluated by minimizing 
a profile likelihood function with respect to the observed 
data and various sources of uncertainty, maintaining all 
correlations among systematic uncertainties [34|. Limits 
on the cross section multiplied by the branching frac- 
tion and the theoretical LQ cross section for /3 — 0.5 are 
shown in Fig. [2] The limit on the LQ mass as a func- 
tion of l3 is determined as shown in Fig .[3l and compared 
to the previous DO Q^, CMS [H H: and ATLAS [H 
results. 
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FIG. 2: (color online) Expected and observed upper limits 
calculated at the 95% C.L. on the LQ cross section as a func- 
tion of Mlq for a scalar leptoquark compared with the NLO 
prediction for /? — 0.5. The NLO cross section is shown for dif- 
ferent choices of the renormalization and factorization scales, 
M = Mlq, h = 0.5 x Mlq, and /i = 2 x Mlq- 



In summary, we have searched for scalar leptoquark 
pair production in the egv^q' final state in 5.4 fb~^ of 
integrated luminosity oi pp collisions at -y/s — 1.96 TeV. 
In the absence of a signal, we exclude the production of 
first generation leptoquarks with Mlq < 326 GcV for 
/3 = 0.5 at the 95% C.L. 
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